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■ Abstract 

CN 

This article completes our study of coherent states in the so-called magnetic-solenoid field (a colinear com- 
bination of a constant uniform magnetic field and Aharonov-Bohm solenoid field) presented in JPA 2010 
and 2011. Here we succeeded to prove nontrivial completeness relations for non-relativistic and relativistic 
l/") ' coherent states in such a field. In addition, we solve here the relevant Stieltjes moment problem and present 

CNl . a comparative analysis of our coherent states and the well-known in the case of pure uniform magnetic field 

Malkin-Man'ko coherent states. 

^ 1 Introduction 
^— > 

A splitting of Landau levels in a superposition of the Aharonov-Bohm (AB) field and a parallel uniform magnetic 
field gives an example of the AB effect for bound states. In what follows, we call such a superposition the 
magnetic-solenoid field (MSF), more precisely MSF is a collinear combination of a constant uniform magnetic 
field of strength B and the AB field, i.e., the field of an infinitely long and infmitesimally thin solenoid with 
a finite constant magnetic flux $. Setting the z-axis along the AB solenoid, the MSF strength takes the form 
B = (0,0, B z ), where 

G\ ' 

OO . $ 

O , B z = B + $(5(x) 5{y) = B H 5(r), B = const, $ = const. (1) 

We use the following electromagnetic potentials!]] A^, assigned to MSF JT]): A = A 3 = 0, and 

o 



/ $ B\ A ( $ B \ 

A * = -^2^ + 2 < A * =X [2^ + 2^ (2) 



with x — rcosO and y = rsmO. Henceforth, for our convenience, we will denote the flux $ as $ = $0(^0 + A*)> 
rS \ where < fi < 1 and <&o = 2TTch/e is the Dirac's fundamental unit of magnetic flux. 

Solutions of the Schrodinger equation with MSF were first studied in [T]. Solutions of relativistic wave 
equations (Klein-Gordon and Dirac ones) with MSF were obtained in [2] and then used in [3] to study AB 
effect in cyclotron and synchrotron radiations. A profound study of these solutions and related problems can 
be found in [U [51 [SJ [71 [HI [S] and [TU]. It is important to stress that in contrast to the pure AB field case, 
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where particles interact with the solenoid for a finite short time, moving in MSF the particles interact with 
solenoid permanently. This opens more possibilities to study such an interaction and correspond a number of 
real physical situations. 

Constructing coherent states (CS) for non-relativistic and relativistic particles in the MSF is a nontrivial 
problem, in particular, due to the non-quadratic structure of particle Hamiltonians in this case. For the first 
time, CS in the MSF were constructed both in non-relativistic and relativistic cases in pTJ [12]. However, some 
problems related to the constructed CS remain still open. In particular, the completeness relations for the CS 
were not presented. In the present article we prove these relations for non-relativistic and relativistic CS in 
MSF. In addition, we solve the relevant Stieltjes moment problem and present a comparative analysis of the CS 
in MSF and the well-known in the case of pure uniform magnetic field Malkin-Man'ko coherent states [13]. 



2 Non-relativistic stationary states 

Let us consider a quantum behavior of a non-relativistic spinless particle with the charge q = — e (e > 0) and 
the mass M in the MSF (see Eq. (TT])) in the direction perpendicular to field B (B > 0), i.e. on the xy-plane. 

As is shown in [llj such a behavior is described by two kinds of wave functions „ 2 (t, 0,r), j — 0, 1, namely 

*%UMr) = e^ t ' K $l na (6,p), P= 2 j-, 1=% i = 0,l, (3) 

which are the eigenfunctions of two commuting to each other operators: the Hamiltonian H± = ^jj (P~ + Py) 
and the angular momentum L z — xp y — yp x , where Pk = Pk + ^A^jc, with p^ = —ihdk, k = x, y, and the vector 
potential A is specified by Eq. $Z§. The eigenvalues of H± and L z are given by £ ni — heB/ (Mc)(ni + 1/2) and 
h(l — Iq), I = 0, ±1, . . . , respectively. The existence of two kinds of states ^n},n 2 {t, 0, r ) is connected with the 
presence of AB field (p 7^ 0) and, what follows from that, the breaking translation symmetry in xy-plane. The 
presence of a non-zero flux $ is also visible in two kinds of functions 0ni,ji 2 (#, p)- 

^nlniiO, P) = Me t{l ~ lo)e I n2 , ni (p), m = m, n 2 = m - / - p, l<0, (4) 
tf&WM = Ue t{l ' la ^~ M I ni , n2 (p), ni =m + l + p, ri2 = m, I > 0, 



which are orthogonal set on the xy-plane |111 112] . J\f — y / j/2n is normalization constant with respect to the 
inner product 



(f,9)±= t(e,p)g(6,p)dxdy=- dp d6f{e,p)g{e,p). (5) 
J 7 Jo Jo 

Here, m = 0, 1, . . . , I n ,m{p) is the Laguerre function |14] that are related to the associated Laguerre polynomials 
L^{p) dES as follows 



(p) 



ml 



r(l + m + a) 



1/2 



1 d m 

e - P /2 pa/ 2 LM ^ = ePp - a e -P p ™+«. (6) 

m! dp 



The radial functions I n ,m(p) were taken regularly as r —> when 1 = 0. It corresponds to a most natural self- 
adjoint extension of the differential symmetric operator H±. Considering a regularized case of a finite-radius 
solenoid one can demonstrate that the zero-radius limit yields such an extension, see [5]. 

We know that the set of the functions ^n},n 2 (t, 0, r) is complete due to the self-adjointness of the H±. 
However, it is useful to show explicitly that these functions satisfy the resolution on unity (completeness relation) 
on ;n/-plane. For this purpose, we introduce the retarded S^lx, x') Green function, which we defined as follows 

S Ie \x,x')=Q(At)S(x,x% S(x,x') - l ^ e - 4£ - A * / ^«„ 2 (e,p)^ 2 (^,p'), (7) 

a 

where At = t — t' and a is integers j, I, and m, where / and m are determined in Eq. (]4]). Q(z) is the Heaviside 
step function. Then the unity resolution for states 4>nl,n 2 {'p 1 p) being written with the help of S rot (x, x') has 
the form 

- i S™\x, x')\ At=Q+ =8{x-x')8{y- y') . (8) 



2 



Note that 5 (x — x') S (y — y') = 7 6(9—9') 5(p—p'). It is convenient to introduce an auxiliary function (x, x') 
by using which S(x,x') is represented as 

00 

s( X ,x') = J2 s i j) M> sP(*>*) = * E e ' i£niAt/h ^lnMp)K[{l 2 ( e '>p')> ( 9 ) 

j=0,l l m=0 

where I < for j = and I > for j = 1. 

Now, by employing the states flU and formula 8.976.5 from [13] . we represent \x,x') as: 



Sr'{x,af) = -r-exp 

47T 



i(Z-Z o )A0-i|l(J + M )Ai 



exp{i(p + ^) C ot[fi 7 At/(2Af)]} / Vpp 7 \ 

sin[^7At/(2M)] T{1+ ^ \i sm[h 1 At/(2M)]J ' [ ' 

where AO = 9 — 9' and I v {z) is the modified Bessel function of the first kind. The upper sign in the index of 
i s related to j = and the lower is for j = 1. The representation (|10p matches with the result obtained 
in [6]. Note that S te!t (x, x') is the integral kernel, it can be changed from convenience considerations by changing 
the integration path in the complex plane of p, p'\ e.g. for p = i£, p' = i£' with real positive £, we have 

dpdp'S re %t',6,e',p,p')f(p)g(p') =i 2 ( dZd?S Te %t\9,6',iZ,iOm)9(i$% (H) 



where /(p), <? (p') are arbitrary integrable functions. Considering the limit At = + in Eq. (|10[) . we can use 
the asymptotic formula (8.451.5) from [14] for the Bessel function. Then, going back to the initial variables, we 
obtain 

S\ 3) {x,x') =i2-e^- l ^5{p-p'). (12) 

At=o+ 2tt 

By using the representation 

1 00 

_ e llAe = S(A0), (13) 

/— — 00 

we verify that relation j8j holds, such that the set of the functions ^>n},n 2 (t, 9, r) is really complete. Note that 
the distribution S(x, x') is not defined for At = 0, that is why the time dependent phase in Eq. ([7]) is important. 

3 Non-relativistic CS 



3.1 CS in MSF 

Following the idea of [11] [12], one has to introduce two kinds (j = 0, 1) of instantaneous CS, which are the 
linear combinations of the states <f>n},n 2 (8, p) given by Eqs. (|4]): 



(9 p) = 1 (0 0) 



-' 1 



m=o V r (l + «i)r(l + ra 2 ) 

The CS are labelled by continuous complex parameters zi and 22 ■ Possible values of n\ and n 2 depend on m, i, 
and j according to Eqs . ([4]) . The normalization constants A/} ( | zi | 2 , | Z2 1 2 ) can be calculated from the overlapping 
formula 

($U) $^") ^ = J , (15) 



where 



^ (0) = Qi-^Vi*! 7 :, y/$z%), Kf^Q,^,^), ^( U ' U ) = E© ^+i(2w), (16) 
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for j = f and z k = z' k , k = 1,2. 

Let us remark that <&i^ Z2 (9, p) do not represent a kind of the Gazeau-Klauder coherent stated (GKSC). 

3.2 Completeness relations 

We are going to prove that CS (|14[) form a complete set on xy-plane, that is it allows a unity resolution with 
the measure dVj(zi,Z2) = W^^zi] 2 , \ z 2\ 2 ) dz\ dz\. This statement is equivalent to the relation 

i=o,i 

f^(x,x')= [ dW^ i d^i^M 2 ) e ^ At/a *i^ a?) 



where W^d^ij 2 , \z2\ 2 ) is the positive weight function and S rct (x,x') satisfies condition ([5]). We consider CS 

defined for almost equal times. We include the time dependent phase e ~ lH ± At / h into the definition of the 
distribution F^(x,x') to provide the consistency of the limit At — » + in Eq. (|T7|) . To prove Eq. (fTT|) . we have 
to find the corresponding weight function Wj 1 ( | zi | 2 , \ z 2\ 2 )- 
Let us check the relations 

F^(x,x') =-i"£sl j \x,x>) j =0,1. (18) 

At_n+ * — ' At=0+ 



At=Q+ 



First, we consider the case j = 0, for which F^°\ after using the explicit form of & Z °^ Z2 (9, p) (see Eqs. (|T4"|)). 
takes the form 

F^Ux ->) T T -** At/ * ^Lyf; "> f^-£j a GK n; f ' fc) (i9) 

^m\n\T(l + m -l-p)T{l-n-k-p) ' 



The auxiliary function G(m, n; Z, fc) is chosen as 

G(m,n;l,k) = f d 2 z,d 2 z 2 ^^ffj, z? "z? z 2 m "'" P '4"'^ 

7 TVodzil 2 , |z 2 r) 



x / rf|z 1 | 2 d|z 2 | 2 |z 1 | m +"| 22 r+"^- fc - 2 ^<(|z 1 | 2 ,|z 2 | 2 ) 
Jo 

/•oo 

= S m , n S hk / dutfou ,n t; ,n - , -''Wff(u,i;), (20) 



where Zk = | ^fe |e' iVJfc (A; = 1,2), u — |zi| 2 , v = |z 2 | 2 and Wq(u,v) — ir 2 Wq(u,v)/ATo(u,v) which is an arbitrary 
positive function that provides Eq. (fT8|) . Taking Wq(u,v) — exp(— u — v) and using the representation T(s) = 
J x 11 ^ 1 e~ x dx of the gamma function, we get 

G(m, n; Z, k) = <5 m ,„ 6 ltk L(l + m) T(l + m - I - fi). (21) 



This function being inserted in Eq. (|19|) gives a correct result for ([18)) with j = 0. 

In the same manner, one can verify the case j = 1. Taking into account ((5J, we see that the validity of f| 1 T[) 
is just the proof of the completeness^of CS. 

We point out that the choice of Wq(u, v) — exp(— u — v) in Eq. ([20)) produces two Stieltjes moment problems 
J °° dx x n W{x) = g{n) = F(l + n), where x and n are respectively taken as u, v and m,m — I — fx. According to 
Pakes's criterion |18) the appeared here Stieltjes moment problems have a unique positive solution e~ x , which 



We recall that GKCS are constracted on the base of a complete set of quantum states in a specific manner, see |16II17| . 
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leads to unambiguous, first time in the literature given, weight function W^ M (|zi| 2 , \z 2 \ 2 ) and at the same to 
unambiguous positive measure di/j(zi, Z\) = W^(|,2i| , \z 2 \ 2 ^)d 2 Z\ d 2 z 2 - 
The weight functions W^{u, v) have the form 

W >, v) = n- 2 e-^ Qi-^CVS, vM, WT(«, w) = Tr"^-^ Q M (V5J, y^)- (22) 

It turns out that W^{u, v) can be expressed via special functions only for /x = and 1/2. The case of p = 
which corresponds to the absence of the AB filed, will be discussed in Section 3.3. In the case p = 1/2, the 
weight functions are 

W] ,2 {u, v) = ^ [erf (y/u + \/v) T erf - v^)] , (23) 

where '— ' is for j = and '+' for j = 1. erf(z) is the "error function" encountered in integrating the normal 
distribution [14] . 

3.3 Zero magnetic flux limit 

Let us study the limit $ = that corresponds to the pure magnetic field without the AB solenoid. 

First of all, we consider such a limit for the stationary states. All topological effects connected with the 
translation symmetry breaking vanish for p = and, in particular, for $ = (lo = 0). As a consequence, the 

shift of the Landau levels is absent for 4>ni,n 2 P) an d it is natural to consider a superposition of j = 1 and 
j = states, 

Cn-i (0, P) + P) = 4>m,l (M, 1 = 0, ±1, ±2, -.. • 

Next, we study the limit of $ = in CS JUJ. Thus, we expect to obtain the Malkin-Man'ko CS [15}. To 
show this, we consider the following superposition of the CS: 

* Z1 , Z2 (0,P) = < /2 *£U(M + < /2 *£ 21 (M- (24) 

At the beginning, we note that the probability distribution of \<fr ZliZ2 (9, p)\ 2 calculated with respect to the inner 
product (f,g)± is equal to 

|* Zl ,, 2 (0,p)| 2 =AAo(|z 1 | 2 ,|z 2 | 2 )+AA 1 (|z 1 | 2 ,|z 2 | 2 ) = e ^l 2 +^l 2 , (25) 

where Afj(\zi\ 2 , \z 2 \ 2 ) = at z^ = z'u are given in Eq. (IT6l) for j = j' . To derive Eq. ((25}l . we employ the 
formula 5.8.3.2 from [T5] and the fact that $ Z {\ Z2 (9, p) arc orthogonal for different j. The density \<& Zl ^ Z2 (9, p)\ 2 

is equal to the normalization constant of Malkin-Man'ko CS, see Eq. (41) in 13 . Then, substituting 3>zi, Z2 {9, p) 
into Eq. (|25]l . we obtain 

m m+|i| oo n r 2 

~ X / TO '( TO +K)' L J nV r l' r 2! 

where r x = m, r 2 = m + \l\. Comparing Eq. (j2"u| with Eq. (41) from [T5], we see that & Zl Z2 (9, p) are just 
Malkin-Man'ko CS. 

Now, let us consider the weight function Wj(u,v) for p — 0. In the limit under consideration, we have 

W°(u,v) + W?(u,v) = tt- 2 , (27) 

where the formula 5.8.3.2 from [TS] was used. Eq. (|27l) is the weight function W°(u, v) for the Malkin-Man'ko 
CS. 
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4 Relativistic stationary states 



Note that relativistic spinless CS are reduced to the non-relativistic case. That is why in the relativistic case, 
only the CS of spinning particles are in a sense nontrivial. In spite of the fact that the algebra of the Dirac 
7-matrices and the spin description in (2 + l)-dim and in (3 + l)-dim are different, considering (3 + l)-dim case, 
we can use technical results obtained for (2 + l)-dim. That is why in the beginning, we consider spinning case 
in (2 + l)-dim. 

The behavior of an electron in MSF in (2 + l)-dim are described by wave functions that obey the Dirac 
equation with such a field, see [5]. These wave functions for given 'polarizations' £ = ±1 (related to one of two 
nonequivalent representation for 7- matrices) and particle/antiparticle energy cpo — ±£± have the form 

^ = e- l{cpat)/h ^J(x\x 2 ). (28) 

In contrast to (3 + l)-dim case, particles and antiparticles in (2 + l)-dim have only one spin polarization states. 
Choosing £ = +1, we deal with 'spin-up' particles, and choosing £ = — 1 with 'spin-down' particles. One can see 
that ipp 1 '(a; 1 ,a; 2 ) = a 2 ^>_2 (a? 1 , a;2 )) where a 2 is a Pauli matrix. That is why we consider here only the case 
£ = 1. The 'spin-up' particle ('+') and antiparticle ('— ') states are denoted as = ^± n, „ ■ The functions 
ip± ni n 2 are common eigenfunctions of the total angular momentum operator J = —ihde + ha 3 /2 and of the 
Hamiltonian H = c(crP± + Mca 3 ). The eigenvalues are equal to H(l — Iq — 1/2) and ±£±, respectively. H 
represent a one-parameter family of self-adjoint Hamiltonians (self-adjoint extensions) that are determined by 
the corresponding boundary conditions. We consider only two special cases: ■& = sign$ = ±1. They correspond 
to a most natural self-adjoint extensions H" 6 . Considering a regularized case of a finite-radius solenoid, one can 
demonstrate that the zero-radius limit yields such extensions, see [5]. The functions ip± n . „ 2 can be represented 
as [12] 

^±L, n2 (M = M jt± , nun2 {a 3 [±n (M) - «rPi] + Mc} u^ ±1 (0,p), 

2 f°° -2 2 - / ~ " 

n (M) = ng— / e - n « T dr, II 2 = M 2 c 2 + (o-Pj. 
V 71 " Jo 

^?,„ 2 , CT (M = €»,.(^)^. Vl = (o)' 

where Pj_ — (Pan Py) an d M.a ± ni n2 are normalization factors with respect to the inner product 

(^VOx^T -1 / dp de^(9,p)^'(9,p). (30) 
Jo Jo 



' 1 = [°X (29) 



For particles the energy spectrum is £ + — [(Mc 2 ) 2 + £±( + i)Y~^ 2 ] for antiparticles the energy spectrum is 

£- = [(Mc 2 ) 2 - 
have the form 



£- = [{Mc 2 ) 2 +£\, _ 1) } 1/2 . The energy £ 2 (a) is given by £ 2 M = 2HceB[m + (l + a)/2\. The functions 4£„ 2 



(+i) J 

(_1)J ' • " KC11C1 6J ly ±(a) la 6 lvcil u y "-i-(ff) 

!]: 

a (6,p) = e^"- 1 ^ 6 I n2 , ni (p), m = m, n 2 =m-l a -p, I < -(1 - 0)/2, 
4!ln 2 A e >P) = eW'-W-* 1 * I nun2 (p), m = m + l a + p, n 2 =m, / > (1 + tf)/2, (31) 

where l a = I — (1 + c)/2, and I m +a,m(p) are given by Eq. ©. These functions form an orthogonal set on 
the semiaxes p > with respect to the scalar product (f,g)±. The existence of two self-adjoint extensions 
is correlated with the irregular behavior of the radial functions I n ,m(p) at the origin when I = and cither 
a = -1 for d = +1 or a = +1 for = -1 0. Note that ft 2 u« n2>±1 = £ 2 ± u^ na>±1 then n (M) ^„ 2 , ±1 = 

£± u nl n 2 ±i- Thus, the spectrum of the operator n (M) is positive defined. We use the n (M) in representation 
(j2"9")l to simplify transition from these stationary states to CS in the next subsection. 

With respect to the self-adjoint operators 1?*, it is know that for any $ the functions i^n},n 2 {0, p) form a 
complete set on the xy-plane. 



G 



To prove this directly and to find an explicit form of the unity resolution, we, similar to the non-relativistic 
case, introduce the retarded Green function S let (x,x') of the Dirac equation, 

S rct (x,x') = 0(At) [S~(x,x') + S + (x,x')] = e(At)[S c (x,x')-S%x,x% 
S c (x,x') = Q(At)S-(x,x')-Q{-At)S + (x,x'), S s = 0(-A*)5 - (a:, x') - Q{At)S+(x, x'), 

S*(x,x>) = .X^ At/ *C^(«,/')C„(«',p'y i (32) 

j,l,m 

where summation in Eq. (|32j) is over all possible quantum numbers j, m, and I is specified in Eq. pip . The 
functions S ,c (x,a;') and S c (x,x') are the causal and anticausal Green functions, respectively. The resolution of 
unity is satisfied if in the following relation holds 



ia 3 S Iet {x,x')\ 



At=0+ 



= 5{x-x') S(y-y') I, 



(33) 



where I is an 2 x 2 identity matrix. We are going to prove that Eq. (|33|) take place in our case. To this end, we 
represent S c (x, x') and S c (x, x') in the form of the Fock-Schwinger proper time integral [6]: 



S c (x,x') 



a 3 (p -(Tp ± ) + Mc A c {x,x'), 



S 5 {x,x') = a 3 (p -<tP_l) + Mc A°(x,x'), 
A c (x 1 x') = / ds f(x,x ! ',s) 3 A c (x 1 x f ) = / ds f(x,x ! ,s), 



(34) 



where Po = an< ^ * ne kernel f(x, x' , s) is given by 



f(x,x',s) = ^2 ^2 fa,i(x,x',s) , f a j(x,x',s) = A a j(s) B a j (s) E 



A<t,i (s) 



<J=±1 i = — oo 

7 



87r 3 / 2 s 1 / 2 sin (7s) 



exp 



i (Mch" 1 ) 2 s + i{l a ~ l Q ) A9 - i {l a + a + p) 7s 



( i (cAi) i . ,, . . 
x exp <^ — — + - (p + p ) cot (7s) 



B„,i (s) = (z) if I + 0, B ff , («) = (2) if 1? = +1, 

B a , (s) = I i±a (z) if 1? = -1, z = e- l7r / 2 V^ 7 /sin( 7 s), S ±1 = (l ± a 3 ) /2 . 



(35) 



The integration path over s is deformed so that it goes slightly below the singular points Sk — kn/^ 1 and — s&, 

k = 1,2, Negative values for s are defined as s — |s|e~ 47T . The kernel f(x,x',s) satisfies the following 

differential equation 



i-f{x,x , ,s) = h- 2 {{McY 



a 3 [pa- erPj 



f(x,x',s). 



(36) 



Remembering that S lct (x,x') is the integral kernel of an integral over the variables p, p' , we can fulfill a 
transformation used above in the non-relativistic case and change the integral path on the complex plane of p, 
p' to a form where p = i£, p' = i£' with real positive £, With respect to irregular behavior of the quantities 
B<jfi (s) as \pp'\ — ► 0, we restrict the range of \pp'\ to < <5 < \pp'\ < 00, with arbitrary S <C 1. Then we take 
the limit s — > + . Under such a condition we use the asymptotic expansion of the Bessel function as s — > + 
and then return to the original variable p, p' on the real semiaxes. Thus, we find that 

lim f a ,t (x, x', s) = !%le« l --W6 (cAt) S(p - p') ~ CT 

for both self-adjoint extensions $ = ±1. Finally, using representation (|13p . we obtain 



lim f(x, x', s) = ichd (cAt) S(x — x) S(y — y)I. 

s-s-0+ 



(37) 
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In the same manner, we can get 



lim f(x,x',s) = ~ichS(cAt) S(x - x') S(y - y')I. 



(38) 



Taking into account that the kernel f(x,x',s) has no any singularity in the lower part of complex plane of 
s, the integral A rct (x, x') = 0(At) [A c (x,x') — A c (x,x')] can be represented as 



A ret (x,a;') = 9(At) / dsf(x,x',s), 



(39) 



where T is a clockwise circle, which connects the points s = +0 and s = • e i7r , and passes in the lower part 
of complex plane of s. If conditions (|37|) and ([3"8| hold, the function 



S ret (x,x')= a 3 (p -<rP_ 



Mc 



A let (x,x') 



(40) 



satisfies Eq. (|3"3"|) and is indeed the retarded Green function of the corresponding Dirac equation [15] . 



5 Relativist ic CS 



Generally speaking, in the relativistic case, the Dirac Hamiltonian is not quadratic in the momenta. Due to this 
fact the time evolution of instantaneous CS on xy-plane (see details in |12j ) is not trivial. However, because the 
time evolution of these states is unitary, it is enough to show that the set of such initial CS is complete. 

For instance, CS for massive spinning ('spin up') particle in MSF on the xy-p\ane and in (2 + 1) dimension 

are 

*8,x,» 3 (M = p [±n (Af) - <xP±] +Mc]n^ Z2!±1 (0,p), u<£, atff (M = *£k<r(Mw<n (41) 



where $ 



(i) 



?, p) are defined in the similar way as the CS of non-relativistic electron, see Eqs. (fT4")l . Taking 



Eq. (|2T)f into account, one can see that the CS (|4"Tj) can be written as 



1 V (6 o) 

m — 71 — m — 12T L^i ±,zi,z 2 \ w > 



E 



' v/rxl + noIXH-na) 



(42) 



where m, n 2 , j and I change according to Eq. (j3"Tj) . The normalization constants A^j,±(|zi| 2 , 1 2:2 1 2 ) can be 
calculated from the overlapping formula 



±,Z1,Z2 ' ±,Z l,z' 2 



2Mc 



D v/^ ± (kil 2 ,k2| 2 )M^ ± (NilM4l 2 ) 



(i) 

Zl,Z 2 ,±l> 



±n (M) + Mc 



*(/) 

z l. z 2 .±! 



where the inner products (-,-) D and (•, •) , are defined by ([30]) and ([5]), respectively. 

Representation (|42l) is like (fl4l) , so that the unity resolution in CS can be done in a form similar to Eq. (flTl) . 
Taking into account the structure of the retarded Green function S Tei (x, x') given by (|32|) . we find 



3=0,1 



At=0+ 



= - ? a 3 S re W) L, n+ , fc = l,2, 



^(1,1')= / d 2 z 1 d 2 z 2 W^(\z 1 \ 2 ,\z 2 \ 2 )J2 e 



-»<cfto(Af)At/R<|v(i) 



*&U(M*$L(*V), (43) 



where the weig ht function Wf(\ Zl \ 2 ,\z 2 \ 2 ) is the same as in (|17[) . As in the non-relativistic case, the consistent 
limit as At — > + can be considered due to the inclusion of an appropriate time-dependent phase factor in the 
definition of F^'(x,x') in (J43j) . In the relativistic case under consideration, we have two such different factors, 
one for particle states, another one for antiparticle states. The proof of the resolution identity (|43]) is quite 
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s(i) 



, p). Using the same weight function W^dzij 2 , | ^2 1 2 ) , we obtain a similar result 



similar to the one for & Zl ,z 2 
for the massless fermions. 

To complete our consideration, we consider (3 + 1) case. The domains of (3 + 1)-Dirac Hamiltonian in MSF 
are trivial extensions of the corresponding domains mentioned in (2 + 1) case, that is why we use for self-adjoint 

(3 + 1)-Dirac Hamiltonian the same notation H® . Of course, in this case = C7 (j2 k=1 2 3 7 fc -P fc + Mc\, 

P 3 — ihd z , and 7 , 7 fe are the 4x4 Dirac gamma matrices. In particular, we consider the CS for spinning 
particle in (3 + l)-dim, which can be constructed by using the set of the orthogonal stationary states defined in 
Eq. (A. 39) from [T2]. The latter can be reduced to ones in (2 + l)-dim as follows 



(3) 

±,s,p 3 ,«i,n 2 



(,:) 



cxp 



--(cll [M)t+p 3 z) 



(i) 

±,S,p3,7li,n 2 



0_l), 



= ±1, 



±,S,P3,m,n2 



M v . 



( M~ x (p 3 /c + sM) + 



^ [M" 1 [f/c+sM) - lj a 3 ^ nun2 (6,p) j 



hp) \ 



(44) 



Here M — \J M 2 + {pz/c) 2 ; p% is z-component of the momentum; flo(M) = IIo(M) „, where fto(M) is 



M=M 



defined in (|2"9"]| ; M P3 is an additional normalization factor; s are eigenvalues of the spin operator S g , 



E Z + S Z H^ )/Mc 2 ; 



and ^ ni , na (e, P ) = ^ nitna (e, P ) 



, where ip± nin2 (@,p) is defined in ([29 



M=M 



Using the matrix structure (|44l) we can find instantaneous CS in (3 + l)-dim in the similar form 



Ci) , 

±,S,P3,Z1,Z2 1 



lp)=M % 



Here the two component column ' z Z2 (9, p) is defined as 



/ 



V 



M 



1 (p 3 /c + sM 
Af-i (p 3 / c + sM 



* 3 *^ ZuZ2 (9,p) 



(45) 



(3) 



(9,p) = * 



±,Zl,Z2 



(0,P) 



M=M 



±,zi,z 2 1 



p) are CS in (2 + l)-dim given by (1421) . Thus, we see that for given p^ and s the unity resolution 



where ^ 

in CS for (3 + l)-dim case is reduced to (2 + l)-dim considered above. Note that representations (f4"4"|) and P5|) 
are convenient for the non-relativistic limit. 

Another types of instantaneous CS in (3 + l)-dim that allow one to construct relativistic time-dependent CS 
were obtained in the work [T^], see Eq. (40) there. They have another matrix structure. In the same work it 
was demonstrated that instantaneous CS in (3 + l)-dim are reduced to ones in (2 + l)-dim. That is why the 
unity resolution in terms of CS in (3 + l)-dim is also reduced to (2 + l)-dim case considered above. This allows 
one to prove that the relativistic time-dependent CS given by Eq. (89) from 12 form a complete system on the 
light cone hypersurface ct — z = const. 
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